Imperial/TP/97-98/8 
Applying the linear S expansion to icfr' 

M. P. Blencowe, H. F. Jones, and A. P. Korte 
The Blackett Laboratory, Imperial College of Science, Technology and Medicine, 
London SW7 2BZ, United Kingdom 
(February 7, 2008) 

Abstract 

The linear 5 expansion (LDE) is applied to the Hamiltonian H = \{p 2 + 
m?x 2 ) + igx 3 , which arises in the study of Lee- Yang zeros in statistical me- 
chanics. Despite being non-Hermitian, this Hamiltonian appears to possess a 
real, positive spectrum. In the LDE, as in perturbation theory, the eigenval- 
ues are naturally real, so a proof of this property devolves on the convergence 
of the expansion. A proof of convergence of a modified version of the LDE is 
provided for the ix 3 potential in zero dimensions. The methods developed in 
zero dimensions are then extended to quantum mechanics, where we provide 

numerical evidence for convergence. 

PACS numbers: 03.65.Ge, ll.15.Tk, ll.10.Jj 



1 



Typeset using REVTpjX 



I. INTRODUCTION 



In a recent paper Bender and Milton carried out an investigation of the following 
non- Hermit ian Hamiltonian: 

H =^(p 2 + m 2 x 2 ) +igx 3 . (1) 

Among other results, they provided analytical evidence for the remarkable property that its 
eigenvalues are real and positive. This Hamiltonian and its field theory counterparts arise 
in the study of the zeros of the Ising model partition function as a function of complex 
magnetic field, the so-called Lee- Yang zeros (see, e.g., Ref. || and references therein). The 
distribution of zeros in the complex parameter plane of the partition function of a given 
system can yield useful information concerning its phase transitions. For the Ising model, 
the zeros lie on the imaginary magnetic field axis and, for imaginary field, the effective theory 
is an interacting scalar field theory with dominant interaction i<j) 3 for small fluctuations || . 

In order to give a proper analytical demonstration of the reality of the eigenvalues of 
Hamiltonian (|]), as well as investigate the higher dimensional quantum field theory ana- 
logues, we require an expansion method which converges for finite, but otherwise arbitrary 
values of the coupling g and mass m. The method must therefore be necessarily nonper- 
turbative in these parameters. Bender and Milton employed a variant of their previous 5 
expansion method [f|, in which the interaction icf) 3 is replaced by (i(j)) 2+s . A Taylor expan- 
sion in 6 of the desired quantity (iV-point Green function, nth. energy eigenvalue, etc.) is 
obtained and then S is set equal to one. The terms of the resulting series are nonperturbative 
in m and g. However, in this method it is extremely difficult to go beyond first order in 5 
and hence test for convergence. 

In the present paper we employ an alternative expansion method, the so-called linear 
5 expansion (LDE). The LDE has been employed as a non-perturbative approximation 
method to study problems in, for example, 4 theory P-fil, quantum chromodynamics [HH, 



relativistic nuclear models (TO] , and electron dynamics in disordered systems |TI[. The 
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method involves constructing a modified action which involves the original action of the 
theory, S, and a soluble trial action, Sq, containing one or more variational parameters Aj: 

S S = S + 5(S-S ). (2) 

The Green function of interest is evaluated as a power series in 5, which is then set equal 
to 1. For the success of the method the trial action Sq needs to be simple, so that one 
can perform high order calculations, but also as close to the true action as possible, so that 
expanding around So is a reasonable procedure with a good chance of converging. Although 
5*5=1 = S is independent of the Aj, there is nonetheless a residual dependence on Aj in the 
truncated series evaluated at 5 = 1, and it is therefore necessary to choose these parameters 
according to some well-defined criterion. Perhaps the most commonly used criterion is the 
principle of minimal sensitivity (PMS), according to which the Aj are chosen to be stationary 
points of the truncated expansion, where the dependence is minimal. Whatever criterion is 
adopted, it is to be applied at each order N in the expansion, so that the Aj become N- 
dependent. This feature of the method is crucial for the convergence of the LDE, which can 
be characterized as an order-dependent split between the bare and the interaction terms in 
the action. In many cases where the fixed split of conventional perturbation theory leads to 
a divergent series, the LDE can be proved to lead to a convergent sequence of approximants. 

In this paper we apply different variants of the LDE to the ix 3 potential in zero dimen- 
sions, where we give a proof of convergence of the expansion for the analogue of the vacuum 
persistence amplitude Z. We then go on to consider the one-dimensional problem and 
provide numerical evidence of convergence of the LDE for the finite-temperature partition 
function and the ground-state energy. 

In Sec. II, which deals with the zero- dimensional problem, we first describe how the con- 
ventional LDE converges, but to the wrong answer. One possible resolution is to split the 
integral up for positive and negative x and apply the PMS separately to each integral |I2" 



Numerically, this gives a sequence of approximants converging to the correct answer, and we 
provide a proof of convergence of this procedure. This LDE variant can be generalized to the 



path integral expression of the 1-D finite-temperature partition function. However, difficul- 
ties are encountered when applied to the ground-state energy, owing to the non-analyticity 
of the splitting procedure. We resolve this problem by using a modified 5 expansion which 
involves a shift parameter. This retains the essential features of the integration splitting 
procedure, but has the advantage of being analytic, thus making the calculations straight- 
forward to carry out and allowing the possibility of generalization to higher dimensions. In 
the last part of Sec. II, a numerical study of the shift method is carried out, and a proof of 
convergence provided using saddle-point techniques. 

In Sec. Ill, where we deal with the quantum mechanical problem, we show that the 
conventional LDE again fails. Numerical evidence of convergence is given for the path 
integral expression of the partition function (using both the splitting and shift techniques) 
and especially for the ground state energy using the shift technique. 

In the conclusion we outline further directions, including proving convergence for the 1-D 
problem and generalizing the shift method to the higher dimensional field theory analogues. 

II. ZERO DIMENSIONS 

The zero-dimensional analogue of the finite-temperature partition function in quantum 
mechanics, or the vacuum persistence amplitude in field theory, is the ordinary integral 

Z = f°° dxe- m2x2+i9x3 '. (3) 

J— oo 

For simplicity we will take m — 0. In that case Z oc g^ 1 ^ 3 and we can set g — 1 without 
loss of generality. In spite of the absence of the convergence factor the integral is still 
well-defined, and can be calculated by splitting up the integration range into x < and 
x > 0, and then rotating the contour by ±7r/6. The exact result, obtained in this way, is 
Z = r(l/3)/ v / 3 = 1.54668588415598, to 15 significant figures. 
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A. Naive application of the LDE 



In evaluating Z using the LDE, the standard approach is to modify the exponent simply 
by adding and subtracting a quadratic term, to give 



The procedure is then to expand Z$ to order 5 , set 5 — 1 and then choose A = A at by 
some criterion or other. In the absence of any additional information the most reasonable 
criterion is the principle of minimal sensitivity (PMS), namely to choose A at as a stationary 
point of Zn- 

The truncated series is J2n=o c m with 



Using (H), Zn can be calculated to high order. However, the residual dependence on A turns 



single maximum, and the value of Z at that maximum steadily approached the exact value 
as N increased. As shown in Fig. 1, where Z^(X) is plotted against A for N = 30, the 
situation here is that there are extremely violent oscillations for small A, which gradually 
decrease in amplitude as A increases, with a final, very broad maximum lying above the 
exact value of Z. 

One's natural choice of a PMS point for A would be the position of this last maximum: 
the residual dependence of Z^ on A around this point is much less than anywhere else on 
the graph. However, it turns out that the sequence Zm(Xn), with Av so chosen, indeed 
converges, but to the wrong answer! The same is true of the previous minimum, and indeed 
none of the stationary points converges to the correct answer. This is a stark warning, in 
this admittedly somewhat pathological case, of the shortcomings of the PMS. It is worth 
mentioning that the scaling behaviour with N of these A at is anomalous: for large N they 
grow like N, rather than the iV 1//3 which would be expected from a saddle-point analysis 




(4) 





out to be radically different from the x A 



case discussed in ||, where, for odd N there was a 
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and which indeed is obtained in the variants of the LDE discussed in the following two 
subsections. In an estimation of the error Z — on the lines of the double saddle-point 
procedure of Ref. |J, the interaction term ix 3 would then be sub-dominant, a clear signal 
that one is on the wrong track. 

B. The splitting procedure 

A possible resolution of this problem is to split up the integration range into positive 
and negative x. This gives Z N (X) = Z^(X) + Z^-(A), where 

Z%{\)= r 'dxe- Xx2 {e Xx2+lx3 } N , (6) 

<j 

and similarly for Z^, where {f{z)}N denotes the Taylor expansion of / to order z . If A is 
taken as A = re - " 1 "/ 3 , where r is real, and the integration contour rotated by x — > xe t7T ^ 6 , 
then Z%{\) = e in/6 Z N (r) where 

Z N (r)= r dxe- rx2 {e rx2 - x3 } N . (7) 

J 

For Z~^ we need to take A = re™/ 3 and rotate the contour in the opposite direction, to obtain 
Z N = e^/ e Z N (r). 

The analysis of is very similar to the x A case discussed in Ref. ||. For odd iV there 
is a single maximum in r, which steadily converges to the exact result. For N — 11 the 
situation is depicted in Fig. 2, where we have multiplied Z N by the appropriate factor of 
^/3. The contrast with Fig. 1 could hardly be more striking. The crucial difference between 
this and the naive application of the LDE is that by splitting the integral we are in fact 
using two different (complex conjugate) values of A^ for the two integrals rather than the 
same value for both. 

The proof of convergence of this procedure is similar to that for the x A case, though 
different in detail. The remainder is essentially given by the integral for the Nth.. 
coefficient in the expansion of Z, namely 
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1 POO 

CN = MJo dxe ~ rx2 ( rx2 - x3 f> (8) 

which can be estimated by saddle-point methods for large N. In this case the PMS value of 
rjy scales as iV 1 / 3 . Writing = aN 1 ^ and x = zN 1 / 3 and using Stirling's approximation 
for N\, the leading behaviour of the integral is 

c N = / dze Nv , (9) 
Jo 

where 

up = -az 2 + ln(a^ 2 - z 3 ) + 1. (10) 

The saddle-point equation dip/dz = is 

2az 3 - 2aV - 3z + 2a = 0, (11) 

which has two positive roots z±. What is required for convergence is that both Ref(z + ) 
and Re(p(z_) be negative, and then the rate of convergence is governed by the larger of 
the two. They become equal at the PMS point, when a is such that Reip(z + ) = Re<p(z-). 
Numerically the solution of this equation is a = 1.0272, when Re(f(z±) = —1.2398. Thus 
the sequence of approximants converges like exp(— 1.23987V). 

C. The shift method 

If we were only concerned with zero dimensions this would be an entirely satisfactory 
resolution of the problem. However, in higher dimensions one is dealing with a functional 
integral, and the generalization of such a splitting procedure is fraught with difficulties, even 
though some progress can be made along these lines, as we show in Sec. IIIB. Fortunately, 
consideration of the quantum mechanical problem suggests another solution, which is im- 
mediately generalizable to higher dimensions, namely to incorporate a linear term, or shift, 
into the (^-modified action. 
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The motivation for introducing such a term is discussed briefly in Sec. IIIC In the 
context of the zero-dimensional model it is worth noting that one of the features of the 
successful stratagem of splitting the integration range and treating each half separately is 
that terms odd in x survive the integration, whereas they cancel in the naive application 
of the ^-expansion. Introducing a linear term into the (^-modified action also avoids such a 
cancellation, but in a simple algebraic way. 

Thus, in zero dimensions the relevant modification of Eq. (f|) is 



where we have introduced a shift ia, which will indeed turn out to be pure imaginary. Notice 
that when 8 = 1 this does reduce to the original integration, as the shift is then immaterial. 

Truncating the expansion at order 8 N and setting 8 = 1 we now have the series Y,n=o c «> 
where 



One can now search for double PMS points in the two parameters a, A. For N > 2 there 
are several such points. For example, at order 16 there are ten solutions, with A lying in the 
range 4-6 and a close to 1. There is no very convincing criterion for choosing between these 
multiple PMS solutions, but they all agree with the exact answer up to the 7th. decimal 
place. 

However, the PMS was only a means to an end, namely to obtain a sequence of ap- 
proximants converging to the exact answer, and in the present case we can instead adopt 
this latter requirement directly as a criterion for choosing a and A. This can be done via a 
saddle-point analysis of the error R N . 

As in the previous subsection, the error is essentially given by the expression for the 
Nth. coefficient cn, which in the present case reads 




(12) 




(13) 




(14) 



S 



The appropriate scaling for a saddle-point approximation in which both parameters play a 
non-trivial role is A = aN 1 ^ 3 , a = fiN 1 / 3 . Then the large- N behaviour of cjv is 



roo 
C N 1 



dze N *, (15) 

J -co 

where now (again setting g = 1) 

ip = -az 2 + \n az 2 + i(z + i(3) 3 +1. (16) 
The saddle-point equation in this case is a quartic: 

2az 4 + 2za(3/3 - a)z 3 - 3(2(3 2 a + l)z 2 - 2i((3 3 a + 3(3 - a)z + 3(3 2 = 0, (17) 

with four complex roots. 

The location of these roots is shown in Fig. |3| for typical values of a and (3, when they 
occur as two pairs of the form ±u + iv. Also shown in that diagram are the paths of 
stationary phase. It is necessary to know the geometry of these paths in order to determine 
which saddle points are encountered, and in what direction, when the contour is distorted 
from its original location along the real axis from — oo to oo. In fact the required path goes 
through all four saddle points A, B, C, D, in each case at a maximum of the integrand 
because of the intervening "sinks" X, Y, Z where the argument of the logarithm vanishes. 

The values of (p at the mirror points A and B are the complex conjugates of those at D 
and C respectively. The criterion for convergence therefore reduces to ensuring that both 
Rep(zc) and Reip(z D ) are negative. Optimal convergence occurs when the two are equal. 

Thus in the parameter space of a and f3 we solve the equation Rep(zc) = Reip(zo) 
and then look for points where this common value is negative. We display the results as 
a contour plot in Fig. |], where the outer contour corresponds to Rep(zc) = and the 
inner one to Rep(zc) = — 1. Any point within the outer contour will give convergence, 
while points within the inner contour will give convergence with an error smaller than e~ N . 
Notice that there is quite a wide range of acceptable values of (3 = aN' 1 ^ 3 , but of course 
this does not include (3 = 0, which would correspond to the naive LDE with no shift term. 
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In Fig. |5| we demonstrate the nature of the numerical convergence up to iV = 30 for 
the optimal values a = 1.79055, (3 = 0.38378. The abolute value of the error is shown, on 
a logarithmic scale. Although this decrease lies within the predicted envelope, the error is 
periodically significantly smaller than the general trend. A similar pattern of convergence 
has been seen ]T3[ in the application of variational perturbation theory to the calculation of 
the strong-coupling coefficients of the quartic oscillator. 

In conclusion, the shift method described in this subsection involves a simple algebraic 
modification of the original integrand and an unambiguous choice of the two variational 
parameters according to the requirement of convergence. The method may be straightfor- 
wardly generalized to higher dimensions, where the calculations involved in the S expansion 
are not significantly more complicated than in ordinary perturbation theory. In the case of 
quantum mechanics the analytic nature of the (^-modified Hamiltonian allows one to go to 
very high orders, as we show in subsection IIIC. 

III. ONE DIMENSION 

We now return to the quantum mechanical problem armed with the insights gained from 
the zero-dimensional analogue. 

The initial evidence for the reality and positivity of the spectrum comes from a numerical 
analysis using the matrix method, whereby the Hamiltonian is rewritten in terms of the 
raising and lowering operators of the bare Hamiltonian and then regarded as a matrix H mn 
in the infinite-dimensional space labelled by the occupation numbers n. If this matrix 
is truncated at some finite order N the resulting spectrum of eigenvalues can easily be 
calculated numerically, up to N of the order of 100. The pattern which emerges from such 
a calculation is that as N increases more and more of the lower eigenvalues emerge as real 
numbers from the amalgamation of complex conjugate pairs. After emergence they each 
tend to a definite positive limit. The ground-state energy is always present from N = 3 
onwards, and for m = g = 1 is stable at E = 0.797342612 to the 9th. decimal place beyond 
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N = 40. 



Since the Ei appear to be real, so also will be the finite-temperature partition function 
= J2i ex v(~/3Ei), which can be calculated to the desired accuracy by including a 
sufficient number of energy levels in the sum. 

To our knowledge there is no proof of the convergence of the matrix method, and our 
ultimate goal will be to provide such a proof for a suitable generalization of the LDE. In 
this section we discuss in turn the three variants of the S expansion which were used in the 
zero-dimensional analogue. 



The standard 5-modification of Eq. ([!]) amounts basically to expanding around a SHO 
of a different, order-dependent frequency according to 



The new effective interaction is igx 3 — Ax 2 , a polynomial. For a polynomial V(x) one has 
the possibility of expanding to very high orders via the use of recursion relations rather than 
the much more cumbersome Rayleigh-Schrodinger perturbation theory. 

These recursion relations are derived by first factoring off the asymptotic behaviour of the 
wave-function according to ip — e~^ y2 x, where y = x^/u and to = (m 2 + 2A)s. Then x an d 
the ground-state energy are both expanded as power series in 5, with x = 1 +Y^=i $ n Xn and 
E = !&>(l+X)J£=i S n e n ). The crucial observation is that with a polynomial interaction of order 
3, the functions Xn are polynomials of order 3n, so that we can write Xn = J2p=i A n)P {iy) p . 
Substituting these expressions into the Schrodinger equation gives coupled recursion relations 
for the e n and A n>p , with e n = 2A„ 2 - 

Using MACSYMA we have evaluated E up to order N = 20. The results are similar to 
those obtained for Z in zero dimensions, namely there are very large oscillations for small 
A, while for larger A there are a number of PMS points, with A scaling like N, where the 
values of E tend to a constant different from the exact value. 



A. Naive application of the LDE 




1 r 



(18) 
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B. The splitting procedure 



In Sec. IIB, we saw how applying LDE separately to Z + and Z~ gave a sequence which 
converged to the correct answer. A natural generalization of the splitting procedure to 
the path integral formulation of the quantum-mechanical partition function is given by 
Z s = Zg + Zg, where 



0/2 



[dx] 6 ± / xdt exp 



(3/2 



13/2 , f 1 
dt < - 

-ft/2 12 



x 2 + (m 2 + 2A±)x 2 + S(igx 3 - \±x 



(19) 



This expression is rather formal as it stands, since it is not possible to evaluate the partition 
function directly in the presence of the step function. A possible resolution [O] is to use the 
integral representation of the step function: 



9{z) 



dq- 



.jiqz 



oo q — le 



(20) 



which leads to 



7± 



2vri 



dq 



Zf{q) 



oo q — it 



(21) 



where 
z t(<l) = I \ dx \ ex P 



ft/2 

dt , 
/2 12 



x 2 + (m 2 + 2A±)x 2 =f iqx + S(igx 3 - \±x 2 ) \ . (22) 



We have checked that the zero-dimensional analogue of ( PT|) gives the same result as splitting 
up the integral directly. 

To first order in 5, Eq. ( [Zip gives (cf. ffM) 

2ig 



z ± = l -z ± h+p 



■(3A± - 1/a; 2 ) 



(23) 



where uj± = m 2 + 2A±, A ± = coth.{\ui±f3) /2ui± ) and Z ± = l/[2 sinh(|a;±/?)] is the partition 
function for the simple harmonic oscillator. 

For each of Z^ there is a single complex PMS point, with A_ = A*]_ and Z~ = (Z + )*. A 
similar calculation can be carried out to order 5 2 , where in each case we find two complex 
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PMS points. In Fig. |] we plot the order 5 and 6 2 approximations to the finite-temperature 
partition function as a function of the coupling g for (3 = 2 and m 2 = 1, and compare 
them with the partition function calculated via Z = J2i e - ^, with the energy eigenvalues^ 
obtained by the matrix method. Both possible solutions at 0(5 2 ) are given, and it will be 
seen that they do not differ appreciably from each other. The graph provides some evidence 
for the convergence of this variant of the LDE method, but it is obviously not conclusive. 
However, it is difficult to go to significantly higher orders in this approach. 



C. The shift method 

While we have seen that the splitting method can be implemented in the calculation of 
Z via the integral representation of the theta function, calculations to higher orders become 
progressively more complicated. Moreover, there is no obvious way to implement the method 
for the calculation of the energy levels themselves. The shift method, in contrast, can be 
applied to both Z and the energy levels, and in the latter case recursion relations allow one 
easily to go to high orders in 5. 

In the context of quantum mechanics, the motivation for introducing a shift arises from 
a modification of the variational approach for the calculation of the ground-state energy 
using a trial Gaussian-like wavefunction 

V> = N[9(x)e- Xix2 + 6{-x)e- Mx2 } (24) 

which distinguishes between x < and x > 0. Because of the symmetry of the Hamiltonian 
under x — > — x combined with complex conjugation it is clear that Ax = A£. 

In the approach of [13[ a delta-modified Hamiltonian is constructed which coincides with 



the Gaussian variational approximation to lowest order in Rayleigh-Schrodinger perturba- 
tion theory, but then allows one systematically to improve on the variational approximation. 



For (3 = 2 sufficient accuracy was attained by including at most 10 energy levels in the sum. 
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The difficulty with this approach in the present problem is the non-analyticity of the Hamil- 
tonian, which considerably complicates the calculation of the higher order terms. 

The shift method can be regarded as a modification of this approach which retains 
its essential features but involves an analytic wave-function and Hamiltonian. Thus with 
Ai = b + ia the wave-function of Eq. Q2"4] ) can be written as ip = ]\f e -( bx2 + mx \ x \) ; with the 
non-analyticity residing in the \x\ factor in the imaginary part of the exponent, which is, as 
a consequence, odd in x. The simplest alternative construction which retains this feature, 
while at the same time being analytic, is ip ~ e -{bx 2 +%ax) _ Thus, a (pure imaginary) shift 
parameter ia has been introduced. A similar shift procedure is in fact described in [[L3[], and 
a method given to ^-modify a Hamiltonian incorporating a shift. Following this method we 
obtain the Hamiltonian 



H& = ^P 2 + Ax 2 - ^mV + ga 3 + 5 



-m 2 (x + ia) 2 — Xx 2 + -m 2 a 2 — ga 3 + ig(x + ia) 3 



(25) 



The corresponding Lagrangian can be used for the calculation of the finite-temperature 
partition function as a functional integral rather than as J2i &~^ Ei - The first-order approxi- 

is 



mation to Z obtained in this way [14 



1 + ?A n ( 3ga + A - ^m 2 



(26) 

where A = coth( v / 2A/?/2)/2 v / 2A and Z = l/[2 sinh( v / 2A/3/2]. We have also calculated the 
second-order approximation. In Fig. [7] we plot both the order 5 and 5 2 approximations to 
Z as a function of the coupling g, where we have again taken m 2 = 1 and (3 = 2. On the 
scale of the diagram the 0(5 2 ) approximation is almost indistinguishable from that obtained 
from the matrix method, which provides strong evidence for the rapid convergence of this 
variant of the LDE for the partition function. 

However, the most convincing numerical evidence for the numerical convergence of the 
shift method comes from a calculation of the ground state energy. This can be done by 
generating recursion relations for Eq. (p5|) in a similar manner to that described in Sec. IIIA. 
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We have applied the recursion relation method to eighth order in 5, taking m = 1 for 
defmiteness. In Fig. |] we show the results of the 0(5) and 0(5 2 ) calculations for a range 
of g, while in Table 1 we illustrate the nature of the numerical convergence up to 0(5 8 ) for 
g — 1. In the latter case we consistently chose the largest PMS value of A. The numerical 
evidence is strong and encourages us to seek a proof of convergence for the one-dimensional 
case. 



N 


A 


a 


Eq 


1 


1.817 


0.3837 


0.7785 


2 


2.169 


0.4445 


0.7957 


3 


2.301 


0.4608 


0.7967 


4 


2.426 


0.4793 


0.79720 


5 


2.525 


0.4930 


0.79729 


6 


2.6107 


0.5040 


0.79732 


7 


2.6865 


0.5135 


0.797334 


8 


2.75475 


0.52145 


0.797339 



Table 1: Results of the shift method for Eq (to be compared with 0.797342612. . . ) 



IV. CONCLUSION 

We have applied the LDE approximation method to the ix 3 potential in zero and in one 
dimension. A proof of convergence was given for the massless case in zero dimensions. The 
methods employed in zero dimensions were then extended to one dimension where numerical 
evidence of convergence for the partition function and the ground state energy was obtained. 
This evidence was particularly strong for the ground-state energy within the shift method. 

The next step is to extend the proof of convergence from zero to one dimension. For Z 
one would need to apply techniques similar to those used in ||, while for the energy levels 



an extension of the methods of Guida et al. [16[ would be necessary. 
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In higher dimensions one has to tackle the problems of regularization and renormalization 
within the LDE. A possible way of doing this has been outlined by the last-quoted authors 
in the context of the calculation of the critical exponents of (p\ using the order-dependent 
mapping method, which is equivalent to the conventional LDE without a shift. It would be 
interesting to generalize their methods to trial actions including a shift, and to apply these 
techniques to the study of Lee- Yang zeros in statistical mechanics. 
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FIGURES 

FIG. 1. Zn for N = 30 as a function of A in the naive LDE in zero dimensions. The dotted 
line shows the exact answer. 

FIG. 2. Z]\f for N = 11 as a function of r using the splitting method in zero dimensions. The 
dotted line shows the exact answer. 

FIG. 3. The location of the saddle points and sinks in the integrand of Eq. (|l5|). The required 
stationary-phase path, passing through all four saddle points, is AX BY C Z D. 

FIG. 4. Contours of Reys in the parameter space of a and j3. Convergence with the shift 
method is assured for all points within the outer contour. For points within the inner contour the 
convergence is faster than e~ N 

FIG. 5. Convergence of the LDE for the shift method in zero dimensions for the optimal values 
of a and (3. The abolute value of the error is plotted on a logarithmic scale against N. 

FIG. 6. The finite temperature partition function versus g, with fixed m 2 = 1 and (3 = 2. The 
solid line is the numerical matrix method approximation, the dashed line the order 5 approximation 
and the dotted lines the order 5 2 approximation using the splitting procedure. 

FIG. 7. The finite temperature partition function versus g, with fixed m 2 = 1 and (3 = 2. The 
solid line is the numerical matrix method approximation, the dashed line the order 5 approximation 
and the dotted line the order 5 2 approximation using the shift procedure. 

FIG. 8. The ground state energy versus g, with fixed m 2 = 1. The solid line is the numerical 
matrix method approximation, the dashed line the order 5 approximation and the two dotted lines 
the order 5 2 approximation with two different PMS values for A, using the shift method. 
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